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Abstract -The  frequency -temperature  (/-T)  behavior  of  a  crystal  res- 
mator  depends  upon  temperature,  and  its  spatial  and  temporal  gradients. 
For  quasi -iso  thermal  changes  the  static  f-T  curve  can  be  used  to  deter¬ 
mine  frequency  shifts  that  occur,  e.g.,  in  oven-controiled  units.  The 
frequency  then  depends  upon  the  parameters  of  the  static  f-T  curve, 
the  temperature  range  over  which  the  oven  cycles,  and  upon  the  oven 
setting  point.  The  maximum  frequency  excursion  has  been  computed 
for  the  A  T  and  SC  cuts  of  quartz  in  terms  of  these  parameters  as  a  func¬ 
tion  of  the  orientation  angle.  When  thermal- transient- compensated 
cuts  are  not  utilized,  oven  cyclings  or  other  temperature  perturbations 
introduce  an  additional  nonnegligible  component  of  the  frequency  shift. 
This  effect  is  quantified  by  means  of  a  simple  mathematical  model. 

The  model  is  capable  of  predicting  the  thermal  transient  effects  for  A  T 
cuts  appearing  in  the  literature.  Simulations,  using  die  model,  disclose 
that  sinusoidal  temperature  variations  with  periods  of  hours  can  readily 
lead  to  frequency  instabilities  much  larger  than  would  be  expected  us¬ 
ing  the  static  f-T  curve  for  the  A  T  cut  This  effect  should  be  greatly 
diminished  in  the  vicinity  of  the  SC  cut. 


Introduction 

IN  THEIR  sixty  years  quartz  resonators  have  shown  them¬ 
selves  capable  of  dramatic  improvements  in  frequency  stabil¬ 
ity.  on  the  average,  one  order  of  magnitude  per  decade.  Be¬ 
cause  of  a  long  history  compared  with  other  solid-state  devices, 
it  is  easy  to  look  upon  this  component  as  the  product  of  a  ma¬ 
ture  technology  and  to  presume  that  the  leveling-off  plateau 
has  been  reached;  that  the  last  bit  of  stability  has  been,  or 
shortly  is  to  be,  wrung  out  of  the  venerable  vibrator. 

Such  a  conclusion  would  be  entirely  wrong.  Prospects  are 
exceedingly  good  that  the  near-term  future  wQl  see  a  further 
hundredfold  improvement  in  the  stability  of  high  precision  os¬ 
cillators,  brought  about  by  using  doubly  rotated  quartz  resona- 
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tors  of  special  orientations.  This  development  fortunately 
coincides  with  stringent  frequency  control  requirements  newly 
imposed  by  the  latest  generations  of  communication  systems 
and  systems  for  navigation  and  position  location. 

In  order  to  identify  an  important  present  contribution  to 
frequency  instability  and  to  separate  out  its  effects,  it  is  con¬ 
venient  to  discuss  first  the  static  frequency-temperature  (f-  T) 
characteristic  of  a  quartz  resonator.  Static  f-T  curves  are  of 
interest  in  their  own  right  in  connection  with  crystals  operat¬ 
ing  without  ovens  over  broad  temperature  ranges,  such  as  in 
temperature  compensated  crystal  oscillators  [2J .  Maximum 
frequency  deviations,  obtained  from  the  static  f-T  curves,  are 
shown  to  be  simply  represented  in  normalized  form  for  such 
applications.  They  are  subsequently  used  for  comparison  with 
the  dynamic  behavior  of  resonators. 

For  crystal  resonators  subjected  to  wide  temperature  varia¬ 
tions,  the  angle  of  cut  is  ordinarily  chosen  to  minimize  the  fre¬ 
quency  excursions  over  the  entire  temperature  range.  High 
precision  oscillators,  on  the  other  hand,  utilize  resonators 
maintained  within  an  oven  system  that  closely  regulates  the 
thermal  environment.  In  this  latter  case  it  is  important  to 
know  how  the  unavoidable  oven  instabilities  affect  the  resona¬ 
tor  frequency  and  how  the  oven  parameters  interact  with  the 
f-T  curve  as  a  function  of  the  crystal  cut. 

Use  of  the  static  f-T  curve  in  conjunction  with  the  known 
characteristics  of  commercial  ovens  leads  to  the  conclusion 
that  frequency  stabilities  orders  of  magnitude  better  than 
those  observed  ought  to  be  possible  with  A  T- cut  resonators. 
The  discrepancy  is  due  to  the  neglect  of  the  dynamic  f-T  he- 
havior  of  the  crystal.  By  making  use  of  various  types  of  dy¬ 
namic  f-T  data  from  the  literature,  a  phenomenological  param¬ 
eter  can  be  extracted  that  satisfactorily  explains  the  effects 
found  in  practice.  When  applied  to  simulations  of  oven  cy¬ 
cling,  it  shows  the  dominance  of  dynamic  effects  in  high  stabil¬ 
ity  situations  when  A  T cuts  are  used. 
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Fig.  1 .  Definition  of  oven  characteristic*  on  static  frequency- 
temperature  resonator  curve.  Cycling  range  is  symmetric  about  set¬ 
ting  point. 


Fig.  2.  Normalized  frequency  excursion  versus  orientation  angle  difference  for  A  T~ cut  quartz  ss  function  of  oven  param¬ 
eters.  Abscissa  reference  angle  is  that  for  which  the  slope  at  the  inflection  point  of  the  static  frequency-temperature 
curve  is  zero.  Angle  varies  somewhat  with  resonator  design  and  construction. 


At  the  SC  cut  the  parameter  linking  dynamic  thermal  behav¬ 
ior  to  frequency  becomes  zero,  and  for  this  doubly  rotated 
quartz  cut  the  f-T  behavior  predicted  from  the  static  curve  is 
capable  of  being  realized. 

In  the  first  part  of  this  paper  the  static  f-T  behavior  of  AT 
and  SC  cuts  is  examined.  The  dynamic  behavior  of  resonators 
is  considered  in  the  second  portion. 

Static  Behavior 

For  most  types  of  quartz  resonators  the  frequency- 
temperature  behavior  takes  the  form  of  the  cubic  curve  de¬ 
picted  in  Fig.  1 .  When  the  resonator  is  operated  in  an  oven 
there  will  be,  in  general,  an  oven  offset  error  A  that  measures 
the  amount  by  which  the  oven  setting  point  misses  the  desired 


reference  point  where  the  first-order  temperature  coefficient 
vanishes.  In  addition,  the  oven  will  have  a  cycling  range  tAT9 
about  the  setting  point,  as  seen  in  Fig.  1 .  As  the  oven  cycles, 
the  frequency  will  vary,  but  in  the  case  shown  the  maximum 
frequency  excursion  will  be  from  the  reference  point  to  the 
frequency  at  the  point  A  +  A 7V  As  A  and  AT0  are  allowed  to 
vary,  the  maximum  frequency  excursion  will  of  course  vary 
and  will  be  a  function  of  A,  AT0,  and  the  shape  of  the  cubic 
curve.  The  cubic  curve  varies  with  the  angle  of  the  cut.  For 
the  <4  f  and  SC  cuts  of  orientations  (YX1)6  =*  35.2  e,(YXwl)+  - 
21 .9 °I0  33.9°,  the  reference  angle  is  taken  as  that  for  which 
turnover  temperatures  coincide  with  the  inflection  tempera¬ 
ture.  In  terms  of  deviations  A8  from  this  angle,  the  maximum 
frequency  excursions  are  shown  for  ^f-cut  quartz  in  Fig.  2. 
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TABLE  I 

Frequency  Deviations  for  A  T-Cvt  Quartz,  with  A8  =  5' and 
Transient  Effects  Omitted 


AO*  S’ 

OVEN  CYCLING  RANGE  (  K  ) 

1 

0  1 

O.OI 

0  001 

OVEN  OFFSET  { K ) 

- 

4.6x10'® 

4.6x10'® 

4. 6x10' 10 

4.6x10'** 

O 

1.4x10"® 

4.7x10'** 

4.7x10'** 

0  01 

1.2x10'® 

1.4x10' 10 

4.7x10'** 

4.7x10'** 

0  001 

1.2x10'* 

1.2x10-*° 

1.4x10' 12 

4.7.10'** 

L 

L 

1.2x10'® 

_ 

1.2x10'*° 

1.2x10'** 

1.2x10'** 

_ 

TABLE  II 

Frequfncy  Deviations  for  A  T-Cut  Quartz,  with  A 8  =  0.5'  and 
Transient  Effects  Omitted 


A  ©  •  .5’ 

OVEN  CYCLING  RANGE  (  K  ) 

_ _ 

1 

0  1 

0.01 

0  001 

*  ( 

1.4.10"® 

1,4.10'*° 

1.4.10'** 

1  “  N  <4elO'9 

*■  'O' 

UJ 

1.6x10*° 

l.'.xlO'1* 

1.6.10  ** 

"  1  _  1 

I  1  °  !i  1.6.13'* 

c  o 

1.6.10  ** 

l  6x10'*® 

100  0 

N  3  A 

1.8x10  ** 

4.6.10*® 

1.6.10  ** 

to 

_ u 

i  sx.c  y  i  a.:o  -1 

1.7.10  *’ 

_ 

5.7.10  1* 

The  angle  difference  is  in  minutes  of  arc,  and  the  ordinate  is 
the  quantity 

|A///|m„/{(Ar0)J  (I  +ipi)J}.  (I) 

Rho  is  defined  as 

p  =  A/ATo  *  2  •  offset/cycle  range.  (2) 

For  Ipl  <  1 ,  a  single  curve  suffices.  The  curves  may  be  approx¬ 
imated  by  the  relation 

O^iK,  ♦*,  Ad)  (Ad)1'5.  (3) 

where  O  is  the  ordinate  expressed  in  ( I )  and  K , ,  K 2  are  simple 
functions  of  Ipl  for  a  given  cut.  Tables  I  and  II  provide  values 
of  | A///|mM  for  A  T- cut  quartz  for  A0  m  5  and  0.S  minutes  of 
arc.  respectively 

The  doubly  rotated  SC  cut  of  orientation  ( YXwi)+  a 
21 .9° /9  2:  ♦33.9°  has  two  angles  to  vary.  Fig.  3  gives  the  vari¬ 
ation  in  O  versus  Ad.  Tables  III  and  IV  present  values  of 
\&flf\mUt  for  Ad  and  A#  *  S',  respectively ,  while  Tablet  V 
and  VI  are  for  Ad  and  A 4  *  0.5',  also  respectively 


Fig.  3.  Normalized  frequency  excursion  versus  orientation  angle  differ¬ 
ence  for  SC- cut  quartz  as  function  of  oven  parameters.  Curves  are 
weak  functions  of  azimuthal  angle  <p. 


TABLE  111 

Frequency  Deviations  for  SC-Cut  Quartz,  with  A0  =  5' 


A®- 

OVEN 

CYCLING  RANGE 

(  K) 

1 

0  1 

0 . 01 

0  001 

_ 

2.1.10'® 

2.1x10'® 

2.1x10-1° 

2.1x10'** 

3d 

H> 

UJ 

o 

6.6.10'® 

2.2x10'*° 

2.2.10'** 

2.2.10'** 

Cft 

It. 

U. 

o 

10  0 

6.4.10  ’ 

t.5.10'** 

2. 2.10' ** 

2.2x10'*® 

z 

Ui 

> 

100  0 

6.4.10'® 

6,6.10'** 

6.6.10'** 

2.2.10'** 

O 

6, 4.10'* 

6.4.10'** 

6.4.10'*® 

6.4.10  ** 

TABLE  IV 

Frequency  Deviations  for  SC-Cut  Quartz,  with  Ad  -  S' 


L0  •  S' 

OVEN 

CYCLING  RANGE  (K) 

_ ! _ 

L  0  1 

0.01 

0  OOi 

X 

- 

6.0x10' ** 

6.0.10'®* 

laJ 

o 

1.8x10'® 

6.5x10'** 

6.5x10'** 

6.5.10  *’ 

Ik. 

Ik 

o 

o 

o 

a 

1.9x10'** 

6.5x10  *> 

6.5x10’®* 

M 

y 

> 

i 

1.6.10'® 

_ 

1.6x10*1 

l.fclO'l® 

6.5x10'®* 

B 

1,6x10'  ** 

1.6x10'®* 

1.6x10'** 

With  reference  to  Table  II,  one  sees  that  an  A  Tout  operated 
in  a  good  oven  with  a  stability  of  a  few  millikelvtns  should 
produce  frequency  stabilities  in  the  order  I0'1*  to  10"M. 

This  in  fact  is  far  from  the  actual  values  obtained. 
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TABLE  V 

Frequency  Deviations  for  SC-Cut  Quartz,  with  A8  =  0.5' 


A9  •  .5* 

OVEN 

CYCLING  RANGE  (K) 

1 

0  1 

O.OI 

ooot 

- 

6. 5x10-5 

6.6x10-5° 

6.Hx10-5I 

6.6x10-52 

*- 

o 

2.0x10-9 

6.7x10-55 

6.7x10-12 

6.7x10-15 

ul 

u. 

Ll 

o 

1.7x10-9 

2.0xl0'u 

6.7x10-55 

6.7x10-59 

o 

o 

z 

III 

o 

o 

1.7x10-9 

1.7x10-55 

2.0x10-55 

6.7x10-55 

> 

o 

o 

1.7x10-9 

1.7x10-55 

1.7x10-55 

1.7x10-55 

TABLE  VI 

Frequency  Deviations  for  SC-Cut  Quartz,  with  A4>  =  0.5' 


O0-.5' 

OVEN 

CrCLINO  RANGE  ( K 1 

1 

0  1 

O.OI 

0  001 

_ 

1.6x10-9 

1.7x10-1° 

1.7x10-55 

1.7x10-52 

* 

Ul 

O 

6.6x10-5° 

2.0x10-55 

2.0x10-'2 

2.0x10-55 

u. 

o 

5.6x10-1° 

6.0x10-12 

2.0x10-15 

2.0x10-5“ 

o 

o 

2 

Ul 

IOO 

5.5x10-*° 

5,1x10-12 

6.1x10-5“ 

Z.OxlO-15 

> 

o 

o 

5.5x10-5° 

5.1x10-52 

5.0x10-5“ 

5.0x10-5° 

Dynamic  Behavior 

That  the  resonance  frequency  of  a  crystal  vibrator  depends 
on  temperature-rate  and  gradient  effects,  as  well  as  tempera¬ 
ture  itself,  has  been  known  for  some  time  [3]  -(26] .  Experi¬ 
mental  evidence  exists  in  a  number  of  forms.  Bistline  [13] 
showed  that  the  f-T cubic  curve  appeared  to  exhibit  cut  angles 
differing  in  AO  by  several  minutes  of  arc  for  different  tempera¬ 
ture  sweep  rates.  Kusters  [20]  showed  how  the  upper  turning 
point  of  an  A  T cut  shifts  in  frequency  with  temperature  rate. 
In  both  cases  the  temperature-time  graphs  were  ramps,  so  that 
t  =  dT/dt  is  a  constant.  The  influence  of  a  “step  function”  of 
temperature  on  A  T cuts  was  shown  by  Warner  [7] ,  ( 1 1  ]  and 
Munn  ( 14) .  For  BT-  cut  and  rota  ted -3T- cut  crystals  [27] ,  the 
frequency  spike  is  reversed  in  sign.  Based  upon  the  observed 
effects,  a  number  of  models  have  been  used  to  characterize 
the  thermal  transient,  or  dynamic  thermal,  behavior  of  quartz 
resonators  (9] ,  (16| ,  (19] ,  (21  ] ,  (24] . 

We  propose  here  a  simple  model  that  explains  the  results  of  a 
variety  of  experiments.  More  careful  experiments,  coupled 
with  a  thorough  analysis  involving  nonlinear  elasticity  theory, 
are  required  to  provide  a  firm  foundation  for  understanding 
dynamic  effects  in  resonators.  It  is  hoped  that  the  results 
presented  here  will  stimulate  such  an  approach  to  the  prob¬ 


lem.  The  traditional  manner  of  expressing  the  static  frequency- 
temperature  behavior  of  quart2  resonators  is  due  to  Bechmann 
[28]: 

A///=  a0  AT  +  b0  AT 2  +  c0  AT3 .  (4) 

In  (4)  AT  =T-  T0,  where  T is  the  present  temperature  and  T0 
is  the  specified  reference  temperature,  often  taken  as  2S°C. 

The  parameters a0,b0,  and  c0  depend  on  material,  cut,  ge¬ 
ometry ,  electroding,  etc.,  but  do  not  depend  on  time,  for  a 
given  resonator  design.  Our  assumption  consists  of  the  inclu¬ 
sion  of  the  term  a  ■  t  in  the  expression  for  the  effective  first- 
order  temperature  coefficient  a(t),  which  now  depends  on  time 
implicitly  through  the  temporal  behavior  of  the  temperature: 

Af(t)lf=  a(t)  ■  AT(t)  +  b0  •  AT(tf  +  c0  •  AT(t)3 ,  (5) 

with 

a(t)  =  a0 +a  ■  t(t).  (6) 

The  quantity  a  is  a  function  of  the  same  entities  as  a0 ,  b0 ,  and 
c0,  as  well  as  several  others  such  as  the  method  of  setting  up 
thermal  gradients  (thickness  gradient,  azimuth-dependent  lat¬ 
eral  gradients),  thermal  conductivity  of  supports  and  elec¬ 
trodes,  etc.,  but  is  a  constant  for  a  given  design.  It  is  treated 
simply  as  a  phenomenological  quantity  to  be  evaluated  sepa¬ 
rately  for  each  crystal  design.  It  is  obvious  that  the  assump¬ 
tion  concerning  <z(f)  could  be  extended  to  the  coefficients  b0 
and  c0  by  introducing  similar  constants  b  and  c;  the  results  are 
qualitatively  the  same,  and  they  will  be  omitted  here. 

Published  Experiments 

For  the  case  of  a  ramp  temperature-time  profile,  where  t  is 
constant,  a(t)  will  reduce  to  a  constant  different  from  o0.  The 
f-T  curve  simply  appears  as  if  it  had  changed  its  apparent  ori¬ 
entation  angle  0.  In  the  case  of  Bistline’s  graphs,  the  value  of  a 
required  to  account  for  the  data  is  +1 .7  X  10~7  s/K1 .  This 
value  is  arrived  at  using  the  known  quantity  da/dO  =  -5.08  X 
10~6/K,  deg  0  for  the  AT  cut  (2) ,  (23] . 

Applying  the  model  to  Kusters’  curves  [20] ,  straightforward 
algebra  yields  the  relation  for  the  normalized  frequency  shift 
of  the  turning  point 

«  =  (rM-  To)- a  t,  (7) 

where  TM  is  the  upper  turnover  temperature  and  T0  is  the  ref¬ 
erence  temperature.  Using  the  parameters  of  the  ^T-cut  crys¬ 
tal  discussed  by  Kusters  [20] ,  the  value  of  a,  for  the  crystal 
used  by  him,  is 

a  -  -2.27  X  10'7  s/KJ.  (8) 

The  Warner  curves  [7] ,  (1 1] ,  (14)  may  be  simulated  closely 
by  making  the  assumption  that  the  thermal  “step  function"  is 
of  the  form  (I  -  e't/T)  and  involves  a  single  system  thermal 
time  constant  r.  The  values  of  a  are  larger  in  magnitude  than 
for  the  Bistline  and  Kusters  experiments.  The  sign  of  a  changes 
according  to  the  electrode  metallization  pattern,  being  nega¬ 
tive  for  thickness  excitation  and  positive  for  electrodes  with  a 
gap  for  lateral  excitation  [11],  It  may  be  expected  that  other 
factors  will  influence  the  sign  and  size  of  i  and  that  a  proper 
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choice  of  resonator  and  electrode  design,  electrode  materials), 
and  method  of  deposition  will  produce  resonators  with  small 
a  values.  The  possibility  of  controlling  a  by  the  use  of  doubly 
rotated  plates  is  discussed  below. 

The  presence  of  a  frequency  spike,  as  found  by  Warner  in  re¬ 
sponse  to  a  thermal  step  function  applied  to  an  A  T-cut  bulk 
wave  resonator,  has  been  recently  verified  for  surface  acoustic 
wave  (SAW)  resonators  of  ST-cut  orientation  by  Parker  [29] . 

Sinusoidal  Temperature  Variation 

Based  upon  the  simple  model  adopted  it  is  possible  to  simu¬ 
late  the  effect  of  sweeping  the  temperature  in  a  sinusoidal 
manner.  It  is  known  that  this  leads  experimentally  to  hystere¬ 
sis  effects.  The  result  of  a  wide  temperature  range  sweep  is 
shown  in  Fig. 4.  The  heavy  line  represents  the  static  f-T curve. 
When  the  temperature  is  swept  as  sin  ut  by  60K  about  the 
30°C  point,  the  resulting  curve  exhibits  an  hysteresis  that  var¬ 
ies  with  the  rate  w  as  shown.  The  effect  is  more  pronounced 
at  the  upper  turning  point  (UTP);  this  is  due  to  the  sign  of  a. 
The  orbits  are  counterclockwise,  also  because  of  the  sign  of  a , 
at  the  UTP. 

The  general  features  of  the  curves  in  Fig.  4,  representing  dy¬ 
namic  thermal  cyclings,  have  been  verified  in  the  main  by  ex¬ 
periments  performed  by  Hughes  [30]  using  sinusoidal  temper¬ 
ature  cyclings.  In  particular,  the  much  greater  deviations  at 
the  UTP  (implying  a  negative  a  in  his  experiment),  compared 
with  the  lower  turning  point,  were  unequivocally  established 
for  the  resonator  examined. 

Fig.  5  shows  the  hysteresis  effects  expected  for  ±10K  excur¬ 
sions  about  the  UTP.  Fora  sweep  with  a  1 0-min  period  the  fre¬ 
quency  change  is  roughly  three  times  that  expected  on  the  ba¬ 
sis  of  the  static  curve.  In  Fig.  6  hysteresis  orbits  for  A  T0  *  5K, 
offset  from  the  UTP  at  by  +5K  and  by  -  10K,  are  shown. 

The  frequency  scale  is  magnified  in  Fig.  7  to  parts  in  10_,° . 
For  temperature  deviations  AT0  of  ±5  X  10-3  K,  the  static 
f-T  curve  is  nearly  a  horizontal  line,  while  the  orbits  for 
the  dynamic  frequency  behavior  are  elliptical,  with  ampli¬ 
tudes  greatly  exceeding  the  static  behavior  when  the  period  of 
the  temperature  cycle  is  shorter  than  about  1  h.  The  constant 
a  from  (8)  has  been  used  in  the  computations.  In  Fig.  8  the 
frequency  scale  has  been  further  magnified  to  the  l(Tn  range, 
with  a  temperature  scale  in  millikelvins  about  the  UTP.  Again 
we  can  see  the  dominance  of  the  dynamic  behavior  compared 
to  the  static,  even  when  the  temperature  variations  take  place 
over  periods  of  the  order  of  a  day. 

Table  VII  provides  values  for  the  fractional  frequency  devia¬ 
tion  |A///|m„  due  to  the  dynamic  thermal  effect  in  /47"-cut 
quartz,  as  a  function  of  both  oven  cycling  range  and  cycling 
time,  assuming  sinusoidal  temperature  variation  and  the  a  de¬ 
termined  from  Kusters’  experiment  [20] .  The  oven  offset  has 
been  taken  as  zero.  The  entries  are  to  be  compared  to  those  in 
Tables  I  and  II.  The  dramatic  influence  of  the  dynamic  com¬ 
ponent  of  the  frequency  shift  can  be  seen  from  the  comparison. 

Gagnepain  [3!  ]  has  confirmed  the  hysteresis  orbits  shown, 
eg.,  in  Figs.  5-8,  for  cyclings  in  a  small  temperature  range  near 
a  turnover  point;  he  uses  a  model  slightly  different  from  that 
of  (5)  and  (6),  where  the  5  term  is  not  multiplied  by  A T(t)  in 
(5);  this  model  has  been  used  by  Anderson  and  Merrill  (9)  and 
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KCKCS  CELSIUS 

Fig.  4.  Simulation  of  frequency  hysteresis  arising  from  sinusoidal  tem¬ 
perature  cycling.  Sixty-kelvin  sweep  about  inflection  temperature. 
Orbital  periods  are  628, 1257,  and  6283  s  for  u  “  10"*,  S  X  10"3, 
and  10,  respectively. 
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SECKES  CELSIUS 

Fig.  5.  Simulation  of  frequency  hysteresis  arising  from  sinusoidal  tem¬ 
perature  cycling.  Ten  kelvin  sweep  about  upper  turnover  point 
Counterclockwise  orbits  stem  from  sign  of  a. 

Koehler  et  al.  [21  ] ,  [24] .  For  small  temperature  sweeps  both 
models  give  equivalent  results,  but  for  wide  sweeps  (5)  predicts 
an  apparent  orientation  angle  shift  for  temperature  ramp  exci¬ 
tation,  whereas  the  modified  model  predicts  a  constant  fre¬ 
quency  offset.  Experimental  results  for  ramp  functions  of 
temperature  have  been  obtained  by  Lagasse  [32] ,  whose  curves 
show  very  interesting  behavior.  Data  were  supplied  on  two 
groups  of  A  T-cut  resonators  having  thicknesses  very  nearly  in 
the  ratio  of  2 : 1 .  For  the  thinner  plates,  60%  showed  only  a 
constant  frequency  offset,  while  the  remainder  showed  both 
an  offset  as  well  as  a  rotation  in  the  apparent  angle.  The 
thicker  plates  showed  the  frequency  offset  alone  in  only  10% 
of  the  cases,  while  90%  had  both  an  offset  and  an  apparent 
angle  change.  Each  resonator  in  a  group  was  of  identical  con¬ 
struction,  but  between  the  groups  the  designs  ware  similar  but 
not  identical. 
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Fjg.  6.  Hysteresii  loops  for  sinusoidsl  temperature  cycles  offset  from 
upper  turnover  point.  Cycling  range  is  five  kelvins;  offsets  are  -10 
and  +5  K.  Orbital  period;  are  2094  s  for  outer  loops  and  6283  s  for 
inner  loops,  corresponding  tow-  3X  10"’  and  1 0-3 ,  respectively. 


Fig.  7.  Elliptical  hysteresis  orbits  for  sinusoidal  temperature  cycling. 
Temperature  sweep  range  S  mQliketvins.  This  figure  is  drawn  for  a  ■ 
-2.27  x  10'7  t/K1.  Orbital  periods  are  3142, 6283,  and  62832  s, 
corresponding  to  u  •  2  X  10'3, 10'3,  and  lO"4,  respectively. 

The  simulation*  seen  in  Figs.  4-8  are  for  sinusoidal  variations 
in  temperature.  In  practice,  the  fluctuations  in  the  thermal 
field  will  have  a  certain  spectral  distribution.  The  result  of  the 
coupling  term  a  will  be  to  produce  a  corresponding  frequency 
spectral  distribution,  with  a  resulting  magnified  frequency 
instability  when  compared  to  what  would  be  expected  in  the 
static  case. 

What  are  the  prospects  for  reducing  a?  Fortunately,  it  ap¬ 
pears  that  by  design  even  xdT-cut  resonators  can  be  produced 
with  smaller  values  of  i.  However,  by  using  doubly  rotated 
crystals  cut  to  orientations  in  the  vicinity  of  the  SC  cut  the 
value  of  a  can  be  reduced  to  a  negligible  value.  These  cuts  are 
compensated  for  a  variety  of  other  effects  as  well  [23] ,  [33] , 
[34] .  For  rotated-Af-cuts  on  the  zero  temperature  coefficient 


Fig.  8.  Elliptical  hyiterem  orbits  for  sinusoidal  temperature  cycling. 
Temperature  sweep  ranges  0.5  and  1  miDikelvin.  Orbital  periods  are 
6.28  X  104  and  6.28  X  10s  s,  for  w  ■  10"4  and  10~5,  respectively.  It 
is  seen  that  even  diurnal  rhythms  contribute  significantly  to  frequency 
instability. 


TABLE  VII 

Frequency  Deviations  due  to  Dynamic  Thermai  Behavior  oe  A  T-Cur 
Resonator  of  [20] 
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locus  (35)  the  approximate  value  fora  is  +1.8  X  10'7  s/K1  for 
EerNis se’s  design.  For  a  given  vibrator  design,  a  plot  of  -a  ver¬ 
sus  0  along  the  upper  zero  temperature  coefficient  locus  in 
quartz  would  look  very  much  like  Fig.  3  in  [33] ;  the  deriva¬ 
tive  da/30  is  zero  at  the  47"  cut  and  reaches  its  maximum  posi¬ 
tive  value  at  the  rotated -AT- cut. 

The  SC  cut  permits  the  possibility  of  a  revolutionary  change 
in  the  stabilities  of  crystal  oscillators,  particularly  in  the 
intermediate-to4ong-term  regimes  With  the  coefficient  a  re¬ 
duced  to  a  negligible  level,  and  with  state-of-the-art  ovens  in 
the  millikeivin  range,  stabilities  in  the  range  10'u  to  10"14 
ought  to  be  achievable.  A  stability  of  6  X  1(T14  for  1 28-s 
sampling  times  has  recently  been  achieved  with  an  SC-cut  crys¬ 
tal  of  special  design,  utilizing  an  airgap  fixture  and  special 
mountings  [36] ,  [37] .  This  stability  also  depended  upon  a 
special  circuit  configuration  in  addition  to  the  SC-cut  resona¬ 
tor  [38] .  The  ability  to  achieve  medium-term  stabilities  on 
this  order  would  make  doubly  routed  quartz  resonators  com- 
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parable  to  rubidium  atomic  frequency  standards  in  precision 
applications.  Medium  precision  stabilities  will  be  available  for 
fast  warmup  oscillators,  and  they  represent  a  second  important 
application  area  for  this  latest  representative  of  a  long  line  of 
quartz  vibrators. 

Conclusion 

The  ac'ual  frequency  instabilities  observed  in  quartz  crystal 
oscillators  for  medium-  and  high-precision  applications  are 
greater  than  can  be  accounted  for  by  the  static  frequency- 
temperature  characteristic  of  the  vibrator  and  the  noise  charac¬ 
teristics  of  the  active  devices  and  ancillary  circuitry.  One 
mechanism  contributing  to  the  instability  is  the  dynamic  f-T 
coupling.  This  paper  has  used  available  experimental  data  with 
a  simple  mathematical  .■  'odel  to  explore  some  of  the  conse¬ 
quences  of  the  dynamic  effect  on  the  characteristics  of  the 
crystal  resonance  frequency. 

Doubly  rotated  quartz  resonators  that  are  compensated  for 
the  dynamic  thermal  effect  offer  promise  in  the  near-term  fu¬ 
ture  of  making  available  very  small,  light  weight,  and  cheap  os¬ 
cillators  with  stabilities  approaching  those  of  rubidium  fre¬ 
quency  standards. 
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